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We use the Gross-Neveu model in 2 < d < 4 as a simple fermionic example for Weinberg's 
asymptotic safety scenario: despite being perturbatively nonrenormalizable, the model defines an 
interacting quantum field theory being valid to arbitrarily high momentum scales owing to the 
existence of a non-Gaufiian fixed point. Using the functional renormalization group, we study the UV 
behavior of the model in both the purely fermionic as well as a partially bosonized language. We show 
that asymptotic safety is realized at non-Gaufiian fixed points in both formulations, the universal 
critical exponents of which we determine quantitatively. The partially bosonized formulation allows 
to make contact to the large- JVf expansion where the model is known to be renormalizable to all- 
orders. In this limit, the fixed-point action as well as all universal critical exponents can be computed 
analytically. As asymptotic safety has become an important scenario for quantizing gravity, our 
description of a well-understood model is meant to provide for an easily accessible and controllable 
example of modern nonperturbative quantum field theory. 



I. INTRODUCTION 

Renormalizability is often described as a seemingly 
technical cornerstone for the construction of admissible 
models in particle physics. Renormalization fixes physi- 
cal parameters of a model to measured values of observ- 
able quantities. A prime physical meaning of renormal- 
izability is the capability of a model to provide an accu- 
rate description of a physical system over a wide range of 
scales at which measurements can be performed. The set 
of physical parameters, say, couplings or mass parameters 
etc., measured at different scales then define the renor- 
malized trajectory in parameter space (theory space). If 
we demand for a specific model to provide a fundamen- 
tal description of nature, the model must be valid on all 
scales, in particular down to arbitrarily short-distance 
scales. The rcnormalizcd trajectory then must exist on 
all scales without developing singularities. 

The requirement of renormalizability can formally be 
verified and realized in perturbatively renormalizable 
theories in a weak coupling expansion. Here, all free pa- 
rameters of a model can be fixed to physical values and 
the renormalized trajectory can be constructed order-by- 
order in a perturbative expansion. This strategy can suc- 
cessfully be applied to theories that exist at least over a 
wide range of scales, but still suffer from a maximum scale 
of UV extension, such as QED [l[ or the standard model 
Higgs sector 0]. If a theory is asymptotically free, i.e., if 
the point of vanishing coupling (Gaufiian fixed point) is a 
UV attractive fixed point, the perturbative construction 
can even be applied on all scales, as in QCD. 

Renormalizability is by no means bound to a pertur- 
bative construction. Even though reliable nonperturba- 
tive information might be difficult to obtain, the concept 
of renormalizability and the existence of a renormalized 
trajectory on all scales can be formulated rather gen- 
erally within Weinberg's asymptotic safety scenario Q. 
Loosely speaking, asymptotic safety is the generalization 
of asymptotic freedom at the Gaufiian fixed point to the 
case of a non-Gaufiian fixed point. As a fixed point of 



the renormalization group (RG) by construction defines 
a point in parameter space where the system becomes 
scale invariant, RG trajectories that hit the fixed point 
towards the ultraviolet (UV) can be extended to arbi- 
trarily high energy/momentum scales, thereby defining a 
fundamental theory, for reviews see 0, [BJ • 

The asymptotic safety scenario has recently become 
an important ansatz for quantizing gravity. In contrast 
to other approaches, this scenario is based on the stan- 
dard gravitational degrees of freedom and also the quan- 
tization procedure proceeds in a rather standard fashion. 
Here, significant progress was made with the aid of the 
functional RG, formulated in terms of a flow equation ]& 
for the effective average action for the metric field [3]. 
In simple truncations, the RG flow of gravity indeed re- 
veals a non-Gaufiian fixed point || - a necessary pre- 
requisite for asymptotic safety. Most importantly, the 
fixed point has remained stable under extensions of the 
truncation, and its universal properties such as the criti- 
cal exponents, in fact, exhibit a quantitative convergence 
under improvements of the approximations involved 0- 
[Tl| . RG relevant directions in theory space have been 
identified and can be associated with a finite number of 
physical parameters to be fixed by experiment. Taken 
together, this provides for a rather strong evidence that 
a quantized version of Einstein gravity can consistently 
be formulated within the asymptotic safety scenario. 

Still many questions are difficult to answer in the con- 
text of quantum gravity, mainly due to technical and 
computational limitations. For a confirmation of the 
asymptotic safety scenario, contact to other nonpertur- 
bative quantization schemes have to be made in a quanti- 
tative manner; first indications for a possible agreement, 
e.g., with dynamical triangulations have already been ob- 
served [l2j]. 

UV-complete scenarios for the matter sector of the 
standard model built on asymptotic safety have also been 
dev elop ed, for instance, for toy models of the Higgs sec- 
tor |13l4l5j also involving non-linear sigma models |16j , or 
for extra dimensional Yang-Mills theories [l7| and gravity 

mta. 
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As asymptotic safety is inherently linked with field the- 
ories in the nonperturbative domain, it appears highly 
worthwhile to identify and investigate other nontrivial 
examples. The present work is devoted to such a detailed 
study. Our benchmark model is given by the standard 
Gross-Neveu model [l9| in d = 3 (or more generally in 
2 < d < 4) spacetime dimensions. This model is known 
to be pcrturbatively nonrenormalizable as the coupling 
constant carries a negative mass dimension. A naive loop 
expansion leads to a series in terms of diagrams with an 
increasing superficial degree of divergence. Proceeding 
in the standard fashion of perturbative renormalization 
would require infinitely many counterterms and thus in- 
finitely many physical parameters to be fixed by experi- 
ment, implying that the theory has no predictive power at 
all. In fact, it has long been known that this conclusion is 
only an artifact of perturbative quantization. By means 
of a Hubbard-Stratonovich transformation, the fermionic 
theory can be partially bosonized such that an alternative 
expansion in terms of the inverse fermion flavor number 
TVf can conveniently be formulated. The large- Nf expan- 
sion turns out to be renormalizable to all orders rather 
similar to a small coupling expansion in a perturbatively 
renormalizable model [20| . 

Whereas this provides strong indications for the exis- 
tence of an interacting Gross-Neveu model in d = 3, it 
remains an open question as to whether this conclusion 
holds for finite Nf. On the other hand, one may wonder 
whether this conclusion about nonperturbative rcnormal- 
izability is indeed profoundly nontrivial: as the partially 
bosonized version of the Gross-Neveu model is identical 
to a Yukawa model. This seems to suggest that the renor- 
malizability in the large- Nf expansion may simply reflect 
the super-renormalizability of the d = 3 Yukawa model. 
In fact, four-fermi models in d — 4 arc known to be re- 
lated to Yukawa models near the Gaufiian fixed point 
[2l|,[22[. In this work, we wish to emphasize that this is, 
in fact, not the case in 2 < d < 4. As we show below, 
also the bosonized Yukawa formulation is renormalized at 
a non-Gaufiian fixed point within the asymptotic safety 
scenario. Similar observations have been made from a 
more pragmatic viewpoint by studying the scaling prop- 
erties of corresponding lattice models towards the con- 
tinuum limit Hj, [H] . 

This work mainly has a pedagogical character. Our 
analysis is performed in a self-contained fashion within 
the modern formulation of the functional RG also to pro- 
vide guidance to the recent literature on asymptotically 
safe quantum gravity. The fixed point analysis and the 
computation of universal properties is performed explic- 
itly and contact is made to the large- Nf expansion where 
fully analytic results for the fixed-point potential and all 
universal critical exponents can be obtained. 

As fermionic models occur in many circumstances of 
particle physics (effective models of QCD) and many- 
body physics (strongly correlated electron systems), our 
analysis of the microscopic completeness of the Gross- 
Neveu model also provides a lesson for the functional RG 



treatment of such systems. Of course, model studies of 
fermionic systems conventionally aim at long-range phe- 
nomena instead of short-distance behavior. As will be 
detailed below, the non-Gaufiian fixed-point facilitating 
UV asymptotic safety in the Gross-Neveu model at the 
same time serves as a quantum critical point associated 
with a 2nd order quantum phase transition towards a 
phase with broken discrete chiral symmetry. This phase 
transition has been studied already with a variety of tech- 
niques [23T - |28j corresponding critical exponents are, of 
course, equivalent to those which we interpret as proper- 
ties of the UV limit of the asymptotically safe model. 

Our work is organized as follows: we begin with sum- 
marizing the essential details of the model in Sect. [II] 
The basics of the asymptotic safety scenario are summa- 
rized in Sect. IIII1 An RG analysis in the fermionic lan- 
guage is performed in Sect. [IVJ Section [V] contains the 
corresponding study in the partially bosonized formula- 
tion including a mean-field and large- Nf analysis and a 
numerical evaluation of the functional RG equations. 



II. GROSS-NEVEU MODEL 

The Gross-Neveu model describes the quantum field 
theory of Nf flavors of masslcss rclativistic fcrmions in d 
space-time dimensions interacting via a four-fermion in- 
teraction term. It allows to study dynamical chiral sym- 
metry breaking. The Euclidean action in d space-time 
dimensions reads 



,x I i=i 



(1) 



where J = J d d x is a shorthand for the integral over the 
d-dimensional Euclidean space-time. The model depends 
on a single parameter which is the coupling constant g 
with mass dimension 2 — d. 

In this work we restrict ourselves to 2 < d < 4; in 
d = 2, the model is asymptotically free and perturba- 
tively renormalizable, as the Gaufiian fixed point is UV 
attractive, d = 4 will turn out to be a marginal case, 
where the asymptotic safety scenario no longer applies 
for integer Nf and the theory becomes "trivial" , i.e., non- 
interacting in the continuum limit, see below and 
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To be specific, we employ a four-component representa- 
tion for the gamma matrices in the d = 3 case in this 
work, i.e., d 7 = 4 where d 7 denotes the dimension of the 
Dirac algebra. In d — 3, the explicit representation of our 
choice for the 4x4 representation of the Dirac algebra 
can be written as 



7o = t 3 <g) t 3 , 7i = r 3 ® Ti , 72 = r 3 <g> t 2 



(2) 



Here, the {t^I's denote the Pauli matrices which satisfy 
nTj = S^tq +\e ljk r k , with k = 1,2,3 and t = 1 2 
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is a 2 x 2 unit matrix. The gamma matrices satisfy the 
anticommutation relation, i.e., the Dirac algebra 



{7m ,lA = 25 P 



(3) 



where /it, v = 0, 1, 2 and I4 denotes the 4x4 unit matrix. 
Moreover we have two additional 4x4 matrices which 
anticommute with all 7^ and with each other: 



73 



-Tl (g> T 



75 = T 2 <8> T , 



7 3 2 = 75 = 14 (4) 



The matrix 735 = i7375 on the other hand commutes with 
7 M and anticommutes with 73 and 75. The action of the 
Gross-Neveu model is invariant under global \J(Nf) trans- 
formations of the fcrmion fields. This also implies invari- 
ance under U(l) 0JVf transformations, i.e., the associated 
U(l)-charge is conserved in each flavor-sector separately. 
The matrix 735 further realizes a U 35 (l) symmetry in 
each flavor sector: 



(5) 



In addition to these two symmetries, the model is sym- 
metric under discrete Z| = {14,75} chiral transforma- 
tions acting on all flavors simultaneously: 



1p l-> -^75 , ip i-> 75^ 



(6) 



(A similar symmetry transformation involving 73 can be 
understood as a combination of Z| and U 35 (l) transfor- 
mations.) 

In continuous dimensions 2 < d < 4, we assume that 
a suitable analytic continuation for the Dirac structure 
exists, such that traces over the algebraic structure yield 
analytic functions in d and d 7 . The chiral symmetry of 
the model can be associated with a Z2 symmetry for 
the order parameter. As we shall discuss below, the in- 
frared regime of the theory is governed by dynamical chi- 
ral symmetry breaking, provided the only parameter of 
the model, namely g, is adjusted accordingly. 



III. ASYMPTOTIC SAFETY AND RG FLOW 
EQUATION 

For a self-contained presentation, let us briefly summa- 
rize the essentials of Weinberg's asymptotic safety sce- 
nario which is based on the underlying general structure 
of the renormalization group (RG). In the space of pa- 
rameters and couplings gi , the RG provides a vector field 
/3, summarizing the RG /? functions for these couplings 
((3)i = (3 9i (gi, g2, ■ ■ ■ ) = dtgi- As the full content of a 
quantum system can be parameterized in terms of gener- 
ating functionals for correlation functions, we can more 
generally study the RG behavior of a generating func- 
tional. Introducing an IR-regulated effective average ac- 
tion Tfe, the RG flow of this action is determined by the 
Wetterich equation @ 



9 t r fe [$] = ^STr{[r( 2) [$] 



Rkj-^dtRk)}, d t = k 



dk' 



(7) 



(2) 

Here, T k is the second functional derivative with respect 
to the field $, representing a collective field variable for 
all bosonic or fermionic degrees of freedom. The func- 
tion Rk denotes a momentum-dependent regulator that 
suppresses IR modes below a momentum scale k. The so- 
lution to the Wetterich equation provides for an RG tra- 
jectory in the space of all action functionals, also known 
as theory space, interpolating between the bare action S\ 
to be quantized Tk-+A — > S\ and the full quantum effec- 
tive action T = T^o, being the generatingfunctional of 
1PI correlation functions; for reviews, see [j| [29l [30j. 

For a fundamental quantum field theory, the RG tra- 
jectory needs to be extendable over all scales which, in 
particular, requires that the UV cutoff A can be sent to 
infinity. This is, for instance, possible if the trajectory ap- 
proaches a point in theory space which is a fixed point un- 
der the RG transformations, i.e., remains invariant under 
the variation of k. Parameterizing the effective average 
action Tk by a possibly infinite set of generalized dimen- 
sionless couplings gi, the Wetterich equation provides us 
with the flow of these couplings dtgi — ft gi (<?i, #2, ■ ■ • )■ A 
fixed point gi^ is defined by 



Pi(9i,*>92,-- 



Vi. 



(8) 



The fixed point is non-Gaufiian, if at least one fixed-point 
coupling is nonzero gj^ ^ 0. If the RG trajectory hits 
a fixed point in the UV, the UV cutoff can safely be 
taken to infinity and the system approaches a conformally 
invariant state for k — > A — > 00. 

For the theory to be predictive, the number of physi- 
cal parameters required for specifying the RG trajectory 
needs to be finite. Considering the linearized flow in the 
fixed-point regime, 



d t gi = Ri 3 {aj,* -9j) + 



9gj 



(9) 



we encounter the stability matrix Bi 3 , which we diago- 
nalize 



(10) 



in terms of right-eigenvectors V/ , enumerated by the in- 
dex /. The resulting critical exponents Q 1 now provide 
for a classification of physical parameters. The solution 
of the coupling flow in the linearized fixed-point regime 
is given by 



9i 



1 



C T V/ 



e 1 



(11) 



where the integration constants C 1 define the initial 
conditions at a reference scale fco- All eigendirections 
with e 1 < are suppressed towards the IR and thus 
are irrelevant. All relevant directions with Q 1 > in- 
crease towards the IR and thus determine the macro- 
scopic physics. For the marginal directions J = 0, 
higher-order terms in the expansion about the fixed point 



4 



have to be regarded. Hence, the number of relevant and 
marginally-relevant directions determines the total num- 
ber of physical parameters to be fixed. The theory is 
predictive if this number is finite. 

For the flow towards the IR, the linearized fixed-point 
flow Eq. © is generally not sufficient and the full non- 
linear (3 functions have to be taken into account. Even 
the parameterization of the effective action in terms of 
the same degrees of freedom in the UV and IR might 
be inappropriate, for instance, if bound states or con- 
densates appear in the IR. This is precisely the case in 
fermionic models beyond criticality; below, we discuss 
asymptotic safety in the Gross-Neveu model therefore 
from both viewpoints, on the one hand in terms of mi- 
croscopic fermionic degrees of freedom and on the other 
hand in terms of a mixed fermionic-bosonic description. 



IV. FERMIONIC FIXED-POINT STRUCTURE 

We begin with a discussion of the fixed-point struc- 
ture of the Gross-Neveu model as it becomes apparent al- 
ready in a very elementary approximation within a purely 
fermionic description. Let us consider only a point-like 
four-fermion interaction, such that our ansatz for the ef- 
fective action reads 



(12) 



where we allowed for a wave-function rcnormalization 
and both and g arc considered to be scale- 
dependent, i.e., a function of fc. This simple ansatz can 
be viewed derivative expansion of the effective ac- 
tion, with the leading-order defined by Z^ =const. This 
expansion can, in fact, be associated with a potentially 
small expansion parameter in terms of the anomalous 
dimension rj^ = —dtlnZ^. Consequently, a running 
wave-function rcnormalization corresponds to a next-to- 
lcading order derivative expansion. Aside from deriva- 
tives, further fermion channels and higher-order interac- 
tions compatible with the symmetries can be taken into 
account. We come back to the role of such interactions 
below. 

Inserting the ansatz (|12|l into the flow equation (J7J, 
the flow of the dimensionlcss renormalized four-fermion 
coupling g, 



(13) 



is given by 

p g = d t g = (d - 2 + 2rH,)g - ^v d if (0) g 2 , (14) 

where v^ 1 = 2 d+1 ir d / 2 T(d/2) and jfy - <D{g). Here, we 
projected the full flow equation straightforwardly onto 
the pointlike limit of the Gross-Neveu coupling; contri- 
butions from further (possibly fluctuation-induced) in- 
teraction channels as well as dependencies on the Fierz 



basis [31H33I] have been ignored for the sake of simplicity. 
The constant if (0) depends on the choice of the regula- 
tor and parameterizes the regulator scheme dependence 
of the RG flow. For instance, for a linear regulator of the 

form laSsi 



1 



-ftfc = Z^r^ip 2 /k 2 ) , r$(x) 



i e(i-x), (15) 



where p 2 = Pq + ■ ■ ■ + pi 2 ,, we have lf(0) = 2/d. Alterna- 
tively, for a sharp cutoff, we find lf(0) = 1. 

Apart from the Gaufiian fixed point we find a second 
non-trivial fixed point for the coupling g which is implic- 
itly given by 



d — 2 + 2 77^ (g* 
4cLyV d l[(p) 



(16) 



At leading order of our derivative expansion we have 
r)^ = and thus 

_ d(d - 2) («j=3) 3tt 2 

^* 8dryVd 4 

for the linear regulator and 

d - 2 (d=3) 7T 2 

* 4:djVd 2 

for the sharp cutoff. The regulator dependence of the 
fixed-point value exemplifies the non-universality of this 
quantity. Nevertheless, the existence of the fixed point 
is a universal statement, as the regulator-dependent con- 
stant if (0) is a positive number for any regulator. More- 
over, the value of the non-Gaufiian fixed point does not 
depend on Nf. In Fig. [TJ we show a sketch of f3 g = dtg. 
The arrows indicate the direction of the flow towards the 
infrared. The theory becomes trivial (non-interacting) in 
the infrared regime for initial values g\ < g*. Choos- 
ing gx > g*, the four-fermion coupling increases rapidly 
towards the infrared and diverges eventually. The diver- 
gence of g at a finite RG scale actually is an artifact of 
the over-simplistic fermionic truncation, but can be as- 
sociated with the onset of chiral symmetry breaking and 
the formation of a fermion condensate. This becomes 
more obvious in the bosonic formulation, see below. In 
fact, the scale for a given IR observable O is set by the 
scale fc cr at which 1/g — > 0: 



O - kt° 



(17) 



where do is the canonical mass dimension of the observ- 
able O. For <7a < g„, the coupling never diverges but 
approaches zero and thus = 0. For <7a > <7* ; we find 



fccr = A 



IffA ~ 9* 
9A 



where the critical exponent is given by 



e 



dg 



d — 2 + 2r;^,„ — 2g, 



(18) 



(19) 
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d t g , 




\g* 







Figure 1: Sketch of the /3 function of the four-fermion cou- 
pling. For g > g, the infrared regime of the Gross-Neveu 
model is governed by chiral symmetry breaking (seemingly 
diverging fermionic interaction g — > oo. For g < g* the model 
becomes a trivial (non-interacting) theory g — s> in the in- 
frared. 

Here, rj^, ^ denotes the value of the fermionic anomalous 
dimension at the fixed point, 77,^* = ?/,/, (.g = .g*). Re- 
lation (fT5)) determines how a given IR observable scales 
when <7a is varied. As the initial condition g\ being 
larger or smaller than g* distinguishes between two dif- 
ferent phases in the long-range limit, the fixed point <?* 
can be viewed as a quantum critical point which divides 
the model into two physically different regimes. 

For our main subject - asymptotic safety -, we ignore 
the infrared physics in the following and concentrate on 
the UV properties induced by the fixed point. In our sim- 
ple fermionic truncation with only one coupling, the sta- 
bility matrix boils down to a number, already represented 
by the critical exponent of the coupling direction in this 
one-dimensional theory space, see Eq. (fl9|) . At leading- 
order derivative expansion, where 77^, = 0, the critical 
exponent is positive for all d > 2, such that the Gross- 
Neveu coupling corresponds to an RG relevant coupling 
being attracted by the non-Gaufiian fixed point towards 
the UV. In this simple truncation, this suggests that the 
Gross-Neveu model can be renormalized and extended 
as a fundamental theory over all scales on RG trajecto- 
ries that emanate from the non-Gaufiian fixed point. As 
there is only one relevant direction, only one physical pa- 
rameter has to be fixed (say the value of the coupling at 
a UV scale, g\) in order to predict all physical quantities 
in the long-range limit. 

It is instructive, to compare these conclusions within 
the asymptotic safety language with standard perturba- 
tion theory near the Gaufiian fixed point <?*,GauB = 0. 
The corresponding critical exponent is 

n d(3 g (g = g*, Gau fi)) , , 

WcauB = ^ = Z - d, (ZV) 

dg 

in agreement with (minus) the naive power-counting di- 
mension of the coupling. At leading-order derivative ex- 
pansion riip ~ and for d > 2, the critical exponent is 
negative and the fixed point thus infrared attractive. A 
UV limit A —> 00 can only be taken if the RG trajec- 



tory emanates from the fixed point, but then the the- 
ory would be noninteracting on all scales and therefore 
"trivial". Within perturbation theory, the conclusion 
is that the Gross-Neveu model is perturbatively 11011- 
renormalizable. Note that this conclusion remains un- 
changed also if the anomalous dimension is taken into ac- 
count: within perturbation theory, 77^ = 0(<?) (actually, 
accidentally 0(<? 2 )), such that 77^ = at the Gaufiian 
fixed point, implying that standard power-counting can 
only be modified logarithmically. 

Let us conclude with a word of caution on the deriva- 
tive expansion in the fermionic truncation: in this sim- 
ple approximation, the fixed-point seems to exist with 
similar properties in any dimension d > 2, in particular 
also in d = 4 and beyond. This conclusion will change, 
once composite bosonic degrees of freedom are taken into 
account. Fluctuations of the latter which are formed 
by fermionic interactions will remove the fixed point in 
the Gross-Neveu model for d > 4 such that no asymp- 
totic safety scenario appears to exist for d > 4 in this 
model. In the fermionic language, the bosonic degrees of 
freedom correspond to specific nonlocal interactions or 
momentum-structures in the fermionic vertices. These 
are not properly resolved in a derivative expansion. As 
d = 4 is a marginal case, the conclusions for fermionic 
theories in d = 4 may depend on the details of the inter- 
action and the algebraic structure of a given model; for 
instance, an asymptotic safety scenario in a standard- 
model- inspired SU(N C ) xU(l) model has been discussed 
in 0. 

A comparison with the asymptotic safety scenario for 
quantum gravity is also instructive: here, the upper crit- 
ical dimension is d = 2 as in the fermionic models, and 
the non-Gaufiian fixed point exist in simple truncations 
based on derivative expansions in all d > 2 [l(| O • ^ is 
tempting to speculate whether strong-coupling phenom- 
ena such as bound-state formation may destabilize the 
fixed point above another so far unknown critical dimen- 
sion. A similar phenomenon has been observed in extra 
dimensional Yang-Mills theories, where a non-Gaufiian 
fixed point exists in 4 + e dimensions but is nonpertur- 
batively destabilized for e > 0(1) [13 • 



V. PARTIALLY BOSONIZED GROSS-NEVEU 
MODEL 

In this section, we study the (UV) fixed-point struc- 
ture of the Gross-Neveu model by employing a (partially) 
bosonized version of the model. We relate our findings 
to the purely fermionic description, study analytically 
the large- Af limit and show how corrections beyond the 
mean-field approximation can be systematically taken 
into account. Formulations of the Gross-Neveu model 
using partial bosonization are used for many aspects of 
the Gross-Neveu model such as the phase structure at 
zero and finite temperature and density [37l - l39j . 

Spontaneous symmetry-breaking and the formation of 
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a fermion condensate in the Gross-Nevcu model can con- 
veniently be studied by introducing an auxiliary field 
a into the functional integral. Formally, we introduce 
this auxiliary field by multiplying the path-integral by a 
suitable Gaufiian factor. This is known as a Hubbard- 
Stratonovich transformation [10, EJ. The partially 
bosonized (PB) action of the Gross-Neveu model then 
reads 



principle, i.e., the gap equation 



^mV + ji (# + iha) i> \ , (21) 



and the functional integral now includes integration mea- 
sures for both fermionic and bosonic fields. As the auxil- 
iary bosonic field only occurs quadratically in the action, 
it can be integrated out again; on the level of the ac- 
tion, this corresponds to replacing the bosonic field by 
its equation of motion, N[fh 2 a = ihifjip. The action 
then reduces again to the Gross-Neveu action upon idcn- 
tifying g = ^j. From the viewpoint of the Hubbard- 
Stratonovich transformation, the Yukawa coupling h is 
redundant as it can be scaled into the sigma field. Only 
the ratio h 2 /in 2 has a physical meaning. In our formula- 
tion, the Yukawa coupling h is implicitly understood to 
carry mass dimension [h] = (4 — d)/2 in order to deal 
with an auxiliary field with canonical mass dimension 
[a] = (d — 2)/2. Moreover, the a field has been scaled 
such that the first and second term in (f2"Tj) are of the 
same order in N{. Under a discrete chiral transformation, 
see Eq. ([5]), the c-field transforms as a i-> —a. From a 
phenomenological point of view a can be considered as 
a bound state of fermions, a ~ ipift. Thus, the vacuum 
expectation value of a is a proper order parameter for 
chiral symmetry breaking, as it is the case in the purely 
fermionic formulation of the model. 



Sa(x) 



s B M 



(23) 



<T (x) 



As we are interested in the UV properties of the model, 
we assume that the ground state is homogeneous. (In 
d = 2, inhomogeneous condensates have been identified 
in some parts of the phase diagram at finite tempera- 
ture and large values of the chemical potential [38|; the 
status for higher-dimensional fermionic models is subject 
to ongoing work, see e. g. [42|, |43|) Such a homogeneous 
ground state cto =const. is then implicitly given by the 
solution of the following equation: 



Co 



p p 2 + h 2 a% 



(24) 



with J p = J nj\d ■ Apparently this equation has a trivial 
solution, (To = 0. Moreover, non-trivial solutions for <7o 
can be found for suitably adjusted values of the four- 
fermion coupling g = h 2 /rh 2 . Since in d — 3 (d = 4) 
the right-hand side of Eq. (f2"4"]l is linearly (quadratically) 
divergent, we impose a sharp UV cutoff A 3> mt 
For d = 3, we find 



h 2 a 2 . 



2 / tt 2 to 2 

7T I TI 2 " 



(25) 



Thus, the fermions acquire a non-zero mass due to the 
spontaneous breakdown of chiral symmetry, provided we 
choose in 2 /h 2 < 2 A/71- 2 . In terms of the four- fermion 
coupling, we can read off a critical value for the dimen- 
sionless coupling g = Ag = Ah 2 /fh 2 above which the IR 
physics is governed by spontaneous symmetry breaking: 



A. Mean-field analysis 

Before we analyze the partially bosonized version of 
the Gross-Neveu model by means of the functional RG, 
let us start with a simple mean-field study, correspond- 
ing to the large- Af limit, in order to rediscover aspects 
of the fermionic language of the preceding section in 
this standard textbook language. Due to the Hubbard- 
Stratonovich transformation, the N{ fermion flavors enter 
only bilinearly and can be integrated out from the corre- 
sponding functional integral. This yields a purely bosonic 
effective theory for the Gross-Neveu model: 



SbW] = I ^m 2 a 2 - A f Trln [# + iha] 
J x 2 



(22) 



9cr= Y 



(26) 



This critical value g CI can be identified with the sharp- 
cutoff value of the non-trivial fixed point g» which we 
found in our study of the fermionic fixed-point structure. 
The role of the critical value as a fixed point becomes ob- 
vious from the fact that if g = g cr the theory is interacting 
but remains massless and ungapped on all scales. From 
the viewpoint of the partially bosonized theory, we find 
that the Yukawa coupling and the boson mass are not in- 
dependent parameters. For a fixed ratio h 2 /fh 2 > g CI the 
IR physics remains unchanged. Thus, the purely bosonic 
description of the theory in this approximation depends 
only on a single parameter as in the fermionic formula- 
tion. 



where Tr denotes a functional trace. Since a depends 
on the space-time coordinates, the action Sb is highly 
non-local and therefore in general difficult to study. The 
ground state (Tq can be obtained from the variational 



B. RG flow of the partially bosonized theory 

Let us now discuss the fixed-point structure of the par- 
tially bosonized theory. A partially bosonized description 
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of the theory is appealing from a field-theoretical point 
as it also forms the basis for the expansion in 1/Nf for a 
large number of flavors. In addition, it allows us to sys- 
tematically resolve parts of the momentum dependence 
of the vertices by means of a derivative expansion. As 
we shall see below, these two expansion schemes are not 
identical and should therefore not be confused with each 
other. For our study we employ the following ansatz for 
the effective action: 

r[^,V,a] = J ^Z a {d^af +i}{Z^ + iha) if, 

+N f U(a 2 )} , (27) 

where we allow all couplings and wave function rcnormal- 
izations Z a $ to be scale dependent. The kinetic term of 
the boson field adds a new aspect: on the one hand, it 
goes beyond the local approximation of simple mean-field 
theory; in terms of the fermionic language, this kinetic 
term corresponds to a specific momentum dependence in 
the scalar s channel of the four-fermion coupling on the 
other hand. As we shall see below, this term and the 
associated wave function rcnormalization receive contri- 
butions to leading order in the large- Nf approximation 
(in order to simplify the large- Nf counting of orders, the 
purely bosonic sector is multiplied by an Nf factor in 
Eq. (J2TJ) similar to Eq. (|2l"])). The large- Nf flow corre- 
sponds to the choice 

Z^\ k ^ A =l, 8^=0. (28) 

This exemplifies the difference between large-iVf and 
derivative expansion, as with this choice we include next- 
to-leading order corrections in terms of a derivative ex- 
pansion in the bosonic sector but treat the fermionic sec- 
tor in the leading-order approximation. 

From our ansatz (|27[) we see that a nonzero homoge- 
neous expectation value for \ha plays the role of a mass 
term for the fermions. By expanding the effective action 
about the homogeneous background a we anticipate that 
condensation occurs only in the homogeneous channel. 
The bosonized Yukawa model is fixed to the fermionic 
Gross-Neveu model by a suitable choice of initial condi- 
tions for r at the UV scale A. This correspondence is 
established by the choice 

^L_ A «1, ^Ifc-A-* 1 ' UK = \m\. (29) 

Thus, the renormalized boson mass m = fh/y/ Z a at 
the UV cutoff A becomes much larger than A and ren- 
ders the boson propagator essentially momentum inde- 
pendent. This compositeness condition for the bosonic 
formulation can be considered as a locality condition at 
the UV scale for the four-fermion coupling g in the purely 
fermionic formulation of the model. 

Since we are interested in the UV fixed-point structure 
of the partially bosonized (chirally-symmetric) Gross- 
Neveu model, we anticipate that a possible non-Gaufjian 



fixed point occurs in the symmetric regime. Therefore, 
we only need to study the RG flow in the symmetric 
regime with vanishing vacuum expectation value for the 
cr field. We then find the following flow equation for the 
dimensionless effective potential: 

dtu = —du + (d — 2 + rj a )u p — 2d 1 Vdl$ (2h 2 p;r]^) 

+-j-2v d l d (u'+2pu";ri*), (30) 

iVf 

where 

u(p) = k- d U(a) , p = ^Z a k 2 ' d a 2 , (31) 

and the dimensionless renormalized Yukawa coupling is 
given by 

h 2 = Z- x Z^ 2 k d -^h 2 . (32) 

The flow equation for the potential can be solved either 
directly or by an expansion around p = 0, 

oo , 

<P) = E ^ ■ (33) 

In this paper, we shall mainly focus on the latter one 
which allows us to directly project onto the (dimension- 
less) bosonic mass parameter A2 and the higher-order 
bosonic couplings \2 n - 
- 2 

A 2 = , X 2n = Z - n k^ d ~ 2n U^\ a=0 . (34) 

The flow equations for the Yukawa coupling as well 
as the anomalous dimensions 77^ = —dt In Z a and 
= — dt In Z^p read 1 

d t h 2 = (d-A + 2r H} +r) a )h 2 

+ ^-8v d hU{ F 1 B)d (0,\ 2 ; n4 ,, Va ), (35) 

iVf 

= s —^- h TO 4 (°;w> ( 36 ) 

1 o Vd u2 (FB)d rr . , \ /o-a 

Vii> = J^ 8 ~J h m i,2 (0,^2;^,^), (37) 

The threshold functions I and m in Eqs. (f3~Tj) and ([33)) - 
(j3"Tl) depend on the details of the regulator. For practical 
computations, we use an optimized regulator [34J - 436I ] for 
the fermionic fields (|T5"j) and the bosonic fields: 

R k = Z a r{p 2 /k 2 ) , r(x) = (~ - lj 6(1 - x) . (38) 

The corresponding threshold functions are listed in 
App. \M These functions essentially describe the thresh- 
old behavior of regularized 1PI diagrams. 



These flow equations agree with those derived in [13 . 1 2SH in the 
symmetric phase, upon a rescaling of the wave function renor- 
malization Z a by a factor of JVf, cf. Eq. I I27I I. 
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Figure 2: Leading-order RG flow in the 1/Nf expansion for the 
partially bosonized Gross-Neveu model in the (h 2 ,\2) plane. 
The red line denotes the critical manifold of points drawn into 
the fixed point towards the infrared. The blue line depicts 
the critical surface attracting the flow towards the IR. The 
critical surface contains all trajectories that emanate from 
the non-Gaufiian fixed point in the UV; its dimensionality 
equals the number of relevant directions and thus the number 
of physical parameters to be fixed. The red dots denote the 
Yukawa Gaufiian fixed point and the non-Gaufiian fixed point, 
respectively. The arrows indicate the direction of flow towards 
the infrared. 



C. RG flow at large N { 

As the nonperturbative renormalizability of the Gross- 
Neveu model beyond d = 2 has been proved to all orders 
in the 1/Nf expansion, it is worthwhile to study asymp- 
totic safety on the basis of the RG flow in the same limit. 
Here, our set of flow equations for the partially bosonized 
Gross-Neveu model reduces to 



d t u 
d t h 2 

Va 



— du + (d — 2 + rj a )u'p 

-2d^v d l[ ) F)d (2h 2 p;r] 4 ,), 
{d-A + 2i li ,+rj a )h 2 1 



^ 2 mf> d (0; %): 



= 0. 



(39) 
(40) 

(41) 
(42) 



dimension is zero. 2 The flow equations for the bosonic 
couplings are essentially driven by the fermion loop, 



d t Xi 



(n(d- 2 + - d)X 2 „ 
-(-l) n n\2 n+2 



(43) 



Fixed point values h 2 , rj a ^, A2,*, A^*, Ag,*, ... can be 
identified as the zeroes of the corresponding j3 functions, 

dth 2 = 0, <9fA2 = 0,... Of course, the Gaufiian fixed 
point with all couplings vanishing solves these fixed point 
equations. As the RG flows for the bosonic couplings 
decouple in the large-iVf limit, a non-trivial fixed point 
requires 7^ 0. This immediately requires 



Va 



(44) 



Whereas this tight relation between the dimensionality 
and the bosonic anomalous dimension here is an arti- 
fact of the large- Nf expansion, a similar relation exists in 
gravity for the graviton anomalous dimension at the fixed 
point as a consequence of background gauge invariance. 
Similar sum rules are known for Yukawa theories with 
chiral symmetries [33| . Such a sum rule for a correspond- 
ing fixed point is also responsible for the universality of 
the BCS-BEC crossover in the broad resonance limit of 
ultracold fermi gases (46J. In the present case, this fixes 
the value of the Yukawa fixed-point coupling, 



(d-4)(d-2) 



d-yVd J (8 — 6d) 



(45) 



The fixed point is interacting for 2 < d < 4 and merges 
with the Gaufiian fixed point in d = 2 and d = 4. Also 
the fixed-point values for the bosonic mass parameter and 
couplings can be given analytically 



A2n,* — 



(-l)"?i!2"+ 2 / d 



(2n — d) \djVd 



n-l 



(d-4)(d-2) 
(8-6d) 



.(46) 



Thus, the fixed-point values for all bosonic vertices A2n 
are non-vanishing. In a purely fermionic formulation of 
the model, these higher bosonic self-interactions corre- 
spond to higher (non-local) fermionic self-interactions. 
In any case, we find that the UV fixed-point theory in 
2 < d < 4 for JVf — > co is not identical to the action Spb 
but involves infinitely many operators. 

It turns out that the alternating series for the full 
renormalized effective potential can be resummed and 
yields a representation of the Gaufiian Hypergeometric 



As the bosonic fluctuations carry no flavor number, we 
observe that 1PI diagrams with at least one inner bosonic 
line decouple completely from the large- Nf RG flow. As 
a consequence, r\ a is non-vanishing in leading order in 
the \jNf expansion, whereas the fermionic anomalous 



2 It is worthwhile to emphasize, that the large- N{ counting is very 
different from scalar O(N) models, where the the anomalous di- 
mensions are zero to leading order at large N{ |44[ . Also the 
structure of the potential equation is very different, such that 
also the search for an exact solution requires a different strategy 
from that of O(N) models [45| . see below. 
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Figure 3: (Resummed) Fixed-point potential ti*(p) in the 
large- Nf limit (blue/solid), see Eq. (|47[) . and low-order poly- 
nomial approximations thereof ranging from p 2 to p 9 for small 
values of p (dashed lines). 



The characteristic polynomial det(_B + 01) of the matrix 
B yielding the eigenvalues — 9 7 via its zeroes is then 
easily found to be 



e-^-d))f[(-Q-(2n-d)) 



(51) 



All large- Nf critical exponents are thus given by d — 4 
and d — 2n. In d = 3 this boils down to one positive 
critical exponent with value 1, i.e., one relevant RG di- 
rection, and infinitely many negative critical exponents 
— 1, — 1, — 3, — 5, — 7, . . . , corresponding to irrelevant RG 
directions. For the case of the Gaufiian Yukawa fixed 
point, the characteristic polynomial of the stability ma- 
trix B is changed to 



-9 - (d - 4)) f[ (-9 - (n(d - 2) - d)) . 



(52) 



function 2-F1 (a, b; c; z) (see e.g. [43]). In terms of renor- 
malizcd fields, the scale invariant fixed-point action then 
describes massless fermions coupled to a scalar boson 
with a potential 



u.(p) = - 



2d 



3d 



;p x 



iFf 



d d (d-4)(d-2) d 

2' ' 2' 6d-8 d^v d ' 



(47) 



This potential has a large-field asymptotic behavior oc 
p? . The small-field region is depicted in Fig. [3] 

The theory has predictive power, as the number of 
physical parameters is determined by the number of RG 
relevant directions corresponding to the number of posi- 
tive critical exponents. At the non-Gaufiian fixed point 
in the large- Nf limit, the stability matrix assumes a par- 
ticularly simple form where only a single column and the 
main diagonal are non- vanishing, 



B 



bfi 2 ,\ 2 b\ 3 , 



h 







2)^2 

b 





A4 , A4 





\ 

) 



(48) 



Therefore only the main diagonal of B enters into the 
stability analysis around the non-Gaufiian fixed point. 
The important non- vanishing entries bh 2 .h 2 an d b\ 2n _\ 2n 
turn out to be completely universal and are given by 



?h 2 .h 2 



d/3 h2 



dh 2 



9=9* 



A-d 



and 



dp* 



d\ 2 



9=9* 



-d + n(d — 2 + 7j CT) *) 



(49) 



2n- d. 
(50) 



As expected, the critical exponents coincide with the 
mass dimension of the Yukawa coupling and the bosonic 
couplings, reproducing simple perturbative power count- 
ing. In total, this yields three relevant RG directions and 
one marginal RG direction. Note that the Gaufiian fixed 
point found in the purely fermionic flow in Sect. lIVl trans- 
lates into a diverging dimensionless renormalized boson 

\!/ 2 -1/2 
mass A 2 ~ g ' . 

Returning to the non-Gaufiian fixed point, we can, of 
course, make contact with the purely fermionic descrip- 
tion and deduce the fixed point of the four-fermion cou- 
pling from the fixed point values of the Yukawa coupling 
and the bosonic mass parameter. We find 



9* 



A 2 * 



d 



d 7 Vd 



d — 2 (d=3) 37r 2 
8 ~' 



(53) 



for the linear regulator which agrees with our findings 
in Sect. IIVI (using the same regulator). In fact, the 
flow equation of the four-fermion interaction g can be 
reconstructed from the flow of h 2 /A2, 



di 



h 2 \ , s (h 2 



8d^v d ( h 2 



-O 



N f 



.(54) 



Using the linear regulator in Eq. (|14p . we observe that 
the flow equation for g and h 2 /A2 are identical in the 
large- Nf limit. Recall that 77^ = in this limit. 

Due to the equivalence of g and h 2 /\2, the quantum 
critical point found in the purely fermionic formulation 
is also present in our study of the partially bosonized 
theory for Nf —> 00, as it should be the case. As it can 
be seen from the scaling law (|18|) . this quantum critical 
point is associated with a vanishing boson mass, i.e., a 
diverging correlation length, in the long-range limit. 

Let us conclude our large- Nf analysis with a word of 
caution on the widely used so-called local potential ap- 
proximation (LPA) in which the running of the wave- 
function rcnormalizations are neglected. If we ignored 
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the running of the wave-function rcnormalization of the 
bosonic field in the present case, the model would artifi- 
cially depend on more than one physical parameter. To 
be more specific, let us consider the mass spectrum of 
the theory in the regime with broken chiral symmetry 
and assume that we have already fixed the mass of the 
fermions. Using the definition of the masses and the flow 
equations of the couplings, we find that the (dimension- 
less) renormalizcd boson mass in the broken regime can 
be written in terms of the (renormalized) fermion mass 
TOf : 



2A 4 po ~ Z'^ippo) 



r7-lT2 2 

Z„ n m 



r i 



(55) 



where po = (1/2)<7q. Neglecting the running of Z a , i. c. 
Z (j = const, as is done in the LPA, we observe that the 
boson mass does not depend on a single physical param- 
eter, as it should be, but on two parameters indepen- 
dently, namely the fermion mass and the (bare) Yukawa 
coupling. By contrast, taking the running of Z a ~ h 2 
into account, the value of the boson mass is fixed solely 
in terms of the fermion mass, in agreement with our fixed- 
point analysis. While this argument might be altered in 
d = 4 space-time dimensions where the Yukawa coupling 
is marginal, it is true for the Gross-Neveu model (as well 
as the Nambu-Jona-Lasinio model) in any dimension d 
in which the flow equation for Z a is non- vanishing even 
at leading order in an expansion in 1/JVf. Therefore the 
flow of Z fj has to be taken into account in a systematic 
and consistent expansion of the flow equations in powers 
of 1/Nf. To be specific, the flow of the order param- 
eter potential (|30p incorporates already fluctuations at 
next-to- leading order in 1 /JVf due to the presence of the 
bosonic loop 3 . However, for a systematic and consistent 
study of the effects of corrections beyond the large-iVf 
expansion the flow of Z a , Z^p ELS well as of the Yukawa 
coupling needs to be taken into account. 



D. RG flow for general flavor number JVf 




Figure 4: Scalar anomalous dimension r\ a for Nf 
2,4,6, ...,100 in d= 3. 



The non-GauBian fixed point of the full Gross-Neveu 
system can now be understood as being sourced from the 
leading large- JVf terms discussed above and the bosonic 
fluctuations inducing a Wilson-Fisher fixed point. De- 
pending on the value of JVf the non-GauBian Gross-Neveu 
fixed point interpolates between the large- Nf fixed point 
for Nf — > oo and the Wilson-Fisher fixed point in the 
formal limit of Nf — > 0. For the latter, the functional RG 
has already proven to be a useful quantitative tool for 
describin g no nperturbative critical phenomena, see e. g. 
Refs. 

Let us repeat the preceding large- JVf analysis, now us- 
ing the full flow equations at next-to-leading order deriva- 
tive expansion, i.e., Eqs. ([3"0|). flSSJ), (|57]l. and ([55]). For all 
quantities of interest, such as critical exponents and fixed 
point values of couplings, a solution of the potential flow 
in a polynomial expansion is sufficient. Confining our nu- 
merical studies to d = 3, all figures are produced within a 
truncated expansion up to 22nd order in cr; quantitative 
results are derived from an expansion to the same order 
in a. In the symmetric regime, a nontrivial fixed point 
in the Yukawa coupling requires the following inequality 
to be satisfied, 



Beyond the limit of large Nf, bosonic fluctuations begin 
to play a role. An immediate consequence is that a new 
fixed point for h = arises for the flow of the effective 
potential for 2 < d < 4. This fixed point of the purely 
bosonic theory is nothing but the Wilson-Fisher fixed 
point which describes critical phenomena in the Ising uni- 
versality class. 



We stress that our parameter Nf plays the role of the number of 
colors iV c in QCD. The number of flavors Nf in, e. g., QCD low- 
energy models is related to the number of involved mesonic scalar 
fields N = 7V f 2 . Thus, a large- N c expansion in QCD models cor- 
responds to a large- JVf expansion in the Gross-Neveu model. For 
an analysis of the role of corrections beyond the large- N c approx- 
imation for the thermodynamics of QCD low-energy models we 
refer to Ref. 0. 



d - 4 + 2t^,* + r? CT) * < 0, for Nf < oo. (56) 

This is because the second term of the Yukawa flow 
Eq. ([SSj) is strictly positive for admissible values of the 
anomalous dimensions r\ a ,r\^ < 0(1)- For instance, in 
d = 3, the sum of the anomalous-dimension terms is al- 
ways slightly smaller than 1, see Figs. 2] and [5j The 
inequality becomes an equality in the large- Nf limit, see 

Eq. dm. 

The resulting fixed point values for the Yukawa cou- 
pling in d = 3 are depicted in Fig. [SJ For increasing 
JVf, the fixed point quickly approaches its large- JVf limit 
(|4"5)) . whereas it tends to zero for small JVf leaving us 
with the pure Wilson-Fisher fixed point of a pure scalar 
model. As the latter is known to exhibit a fixed-point 
potential in the broken regime, i.e., with a non- vanishing 
expectation value of the scalar field cr, we expect such 
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Figure 5: Scalar anomalous dimension for Nf 
2,4,6, ...,100 ind = 3. 




Nf 



Figure 6: Yukawa fixed-point value over Nf = 2, 4, 6, . . . , 100 
in d = 3. 
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Figure 7: Fixed-point potentials it, (p) for Nf = 2, 4, 6, . . . , 100 
in d = 3. The large-iVf result is shown as a thick/blue line. 
Incidentally, for Ni = the potential approaches the Wilson- 
Fisher fixed-point potential with a nontrivial minimum near 
p ~ 0.3 (not shown). 
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a fixed-point potential featuring a nontrivial minimum 
to occur for small Nf. For all integer values of Nf > 1 
Dirac (four-component) fermions, we still observe fixed- 
point potentials in the symmetric regime, in agreement 
with [28j | . Nevertheless, the fixed point seems to occur in 
the broken regime for the model with one two-component 
fermion (corresponding to Nf = 1/2 in our language) (28j |. 

The fixed-point potential in d = 3 for various values 
of Nf is plotted in Fig. [7] in a 22nd-order approximation. 
Also the potential converges rapidly to the large- Nf result 
for increasing values of Nf. 

Let us now turn to the leading universal critical expo- 
nents O 1 ' 2 . The convergence of the polynomial expansion 
of the potential is demonstrated in Tabs. 1 1 1 J 11 where the 
leading critical exponents for Nf = 2 and Nf = 12 are 
listed for increasing truncation order. In each case, the 
leading exponents converge to a stable value. We ob- 
serve a more rapid convergence for larger values of Nf. 
The somewhat slower convergence for more scalar dom- 
inated models is familiar from pure O(N) models. The 
leading two critical exponents are plotted as a function of 
Nf in Figs. |8]and |U Table Hill lists the leading exponents 
for increasing number of Nf, illustrating the approach to 
the analytical large- Nf results. 



Figure 8: Relevant critical exponent O 1 for Nf = 
2, 4, 6, ... , 100 in d — 3. The non-monotonic behavior for 
small Nf is expected, since has to approach O 1 ~ 1.6 
for Nf — corresponding to a correlation length exponent 
v = 1/0 ~ 0.63 of the Ising universality class. 



0- 




N 



Figure 9: First subleading irrelevant critical exponent O for 
Nf = 2, 4, 6, . . . , 100 in d = 3. 
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2n 


e 1 


e 2 


e 3 


e 4 


e 5 


e 6 


4 


0.9928 


-0.8687 


-1.5743 


- 


- 


- 


6 


0.9766 


-0.8743 


-1.0624 


-5.4313 


- 


- 


8 


0.9831 


-0.8721 


-1.0790 


-3.5622 


-10.5959 


- 


10 


0.9821 


-0.8720 


-1.0999 


-3.4194 


-6.8111 


-17.4807 


12 


0.9819 


-0.8723 


-1.0897 


-3.5628 


-5.7949 


-11.2156 


14 


0.9821 


-0.8722 


-1.0911 


-3.5104 


-6.0610 


-8.6408 


16 


0.9820 


-0.8722 


-1.0920 


-3.5062 


-6.1190 


-8.3598 


18 


0.9820 


-0.8722 


-1.0914 


-3.5202 


-5.9849 


-8.9972 


20 


0.9821 


-0.8722 


-1.0915 


-3.5132 


-6.0516 


-8.5869 


22 


0.9821 


-0.8722 


-1.0916 


-3.5135 


-6.0514 


-8.5820 



Table I: Non-Gaufiian critical exponents in d = 3 for increasing polynomial truncations for Nf = 2. The results for the critical 
exponent Q 1 agree within the error bars with the result from Monte-Carlo (MC) simulations [24|, 1/Gmc = Um c ~ 1-00(4). 



2n 


e 1 


e 2 


e 3 


e 4 


e 5 


e 6 


4 


0.9898 


-0.9735 


-1.0701 








6 


0.9903 


-0.9735 


-1.0489 


-3.2583 






8 


0.9903 


-0.9735 


-1.0507 


-3.1714 


-5.5889 




10 


0.9903 


-0.9735 


-1.0505 


-3.1821 


-5.3368 


-8.1011 



Table II: Non-Gaufiian critical exponents in d = 3 for increasing polynomial truncations for Nf = 12. 



N t 


e 1 


e 2 


e 3 


e 4 


e 5 


e 6 


2 


0.9821 


-0.8722 


-1.0916 


-3.5135 


-6.0514 


-8.5820 


4 


0.9775 


-0.9240 


-1.1010 


-3.3910 


-5.7739 


-8.2429 


12 


0.9903 


-0.9735 


-1.0506 


-3.1810 


-5.3665 


-7.6004 


50 


0.9975 


-0.9936 


-1.0143 


-3.0510 


-5.1062 


-7.1789 


100 


0.9987 


-0.9968 


-1.0073 


-3.0263 


-5.0550 


-7.0934 


oo 


1 


-1 


-1 


-3 


-5 


-7 



Table III: Non-Gaufiian critical exponents for various flavor numbers Nf in d = 3 in the p 11 truncation. 



This approach is also visible in the anomalous dimen- 
sion and the fixed point values for the coupling, see 
Tab. IIVI Whereas the fixed-point couplings are non- 
universal (holding for the linear regulator in this case), 
the anomalous dimensions are universal and illustrate the 
inequality (|56p . 

Quantitatively, our results for the leading critical ex- 
ponent can be compared to those studies of the quan- 
tum critical phase transition aiming at the long-range 
physics, as 9 1 is related to the correlation length expo- 
nent v by v = I/O . Together with the scalar anoma- 
lous dimension and corresponding scaling and hyperscal- 
ing relation, all thermodynamic exponents of the phase 
transition are determined. Wherever comparable, our re- 
sults agree quantitatively with the functional RG study 
of [28j where both a polynomial expansion as well as a 
grid solution of the potential was used (note that our Nf 
counts four-component fermions, whereas [28| uses two- 
component fermions). Also the agreement with results 
from other methods such as 1/iVf expansions (2(| and 
Monte Carlo simulations [HI, [24[ is satisfactory. Discrep- 
ancies are mainly visible only in the anomalous dimen- 



sions for small Nf, a feature familiar from scalar mod- 
els. Also the polynomial expansion of the potential con- 
verges somewhat slower for the special case Nf = 1 , where 
our results for the leading exponents are compatible with 
those of 28[ , but subleading exponents seem to require a 
precise solution of the potential flow. 



To summarize, the nonperturbative features of the 
Gross-Neveu model near criticality can well be described 
by the functional renormalization group, as the model in- 
terpolates between two well-accessible limits: the large- 
Nf limit for Nf — > oo and the Wilson-Fisher fixed point in 
the Ising universality class for Nf — s- 0. Our results sug- 
gest that the Gross-Neveu model is asymptotically safe 
for all Nf > and that the model depends on only one 
parameter, even when we take into account corrections 
beyond the large-TVf limit. This might provide helpful 
information, e.g., for a systematic study of the finite- 
temperature phase diagram of the Gross-Neveu model 
beyond the large- Nf approximation. 



13 



Nt 


r?<r,* 




hi 


A 2l * 


2 


0.7596 


0.0320 


4.5565 


0.3956 


4 


0.8870 


0.0138 


5.2820 


0.5846 


12 


0.9644 


0.0041 


5.7206 


0.7263 


50 


0.9917 


0.0009 


5.8746 


0.7822 


100 


0.9958 


0.0005 


5.8983 


0.7911 


oo 


1 





5.9218 


0.8 



Table IV: Non-Gaufiian fixed-point values of the universal anomalous dimensions for various flavor numbers Nf in d — 3. The 
(non-universal) fixed-point couplings hold for the linear regulator. In Monte-Carlo (MC) simulations [24| rff^ — 0.754(8) has 
been found for Nf — 4 two-component fermions (corresponding to Nf — 2 in our language). 



VI. CONCLUSIONS AND OUTLOOK 

We have used the functional RG to describe the UV 
behavior of the Gross-Neveu model in 2 < d < 4 di- 
mensions. In agreement with many earlier results in the 
literature, the model is nonpcrturbatively renormalizable 
by means of a non-Gaufiian fixed point, providing a sim- 
ple example of asymptotic safety in Weinberg's sense. 
The perturbative conclusion about nonrenormalizability 
is a mere artifact of naive power-counting which is only 
justified near the Gaufiian fixed point. 

In this work, we have summarized these conclusions in 
the functional RG language as it is also extensively used 
recently to explore the possibility of quantizing gravity 
within pure quantum field theory. This pedagogic charac- 
ter of our work is amended by new results in the large- Nf 
limit, where the fixed-point potential as well as all critical 
exponents can be computed analytically. Moreover, we 
have also provided finite Nf results, demonstrating that 
the model interpolates between the large Nf limit and a 
purely bosonic model in the Ising universality class. 

Some final comments arc in order: the simple fcrmionic 
Gross-Neveu action is minimalistic in the sense that it 
suffices to put the system into the right "Gross-Neveu 
universality class" . Comparing the fermionic action with 
the fixed-point action in the partially bosonized version, 
we have to conclude that the fixed-point action in the 
fermionic language is far more complicated than the sim- 
ple ansatz (|12|) . In general, it will contain higher-order 
as well as non-local interaction terms. The deviations 
from the simple Gross-Neveu structure, however, are ir- 
relevant operators which do not modify the predictive 
power of the Gross-Neveu model. 

The simple Gross-Neveu action is also incomplete in 
the sense that it does not exhibit all possible four-fermi 
terms compatible with the defining symmetries of the 
model. For instance, a Thirring-like interaction (ipj^ip) 2 
is also invariant under the symmetries of the Gross-Neveu 
model and will thus generically be generated by the RG 
flow. A more complete RG analysis thus has to include 
these terms facilitating the appearance of further fixed- 
points, as is known for the Thirring model (HJ. If so, 
this implies the possible existence of further UV comple- 
tions of fermionic models with the same symmetries as 



the Gross-Neveu model. 

Finally, let as mention once more that the property 
of asymptotic safety in the Gross-Neveu model is tightly 
related to the occurrence of a quantum phase transition 
of 2nd order separating a disordered phase from a phase 
with broken (discrete) chiral symmetry. More generally, 
models with such 2nd-ordcr quantum phase transitions 
are guaranteed to be asymptotically safe. Whether the 
converse is true, i.e., whether asymptotically safe models 
always exhibit a physically relevant order-disorder quan- 
tum phase transition, is an interesting question for future 
studies. 
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Appendix A: Threshold functions 

The regulator dependence of the flow equations is 
encoded in dimcnsionlcss threshold functions which 
arise from 1PI diagrams incorporating bosonic and/or 
fermionic fields. In this work we have employed a linear 
regulator dHH, see Eqs. $T5§ and (ggj. For the thresh- 
old functions appearing in the flow equations in Sect. fVl 
we then find 



AF)d ( 2 ( Tfy \ 1 



,(FB)d, 

'1,1 (Wl,W2;i7V,7?oO = 3 



d(l + wi)(l+w 2 ) 



x 1- 



d + ij l + wi 



l - 



d+2J l + w 2 J ' 
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I _ „ m \2 independent of r]^. 



(F)d, s 1 . 1 ~ Vip 



(1 + cj) 4 d-2 (l + w) s 
•1-^ l\ 1 



2d -4 4/ (1+w) 2 ' 



{FB)d 1 \ 

»7<7 



d+iy (i + w 1 )(i + w 2 ) : 



Note that the linear regulators (|T5|) and (|38|) render 
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